(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 13.2' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     79222,       2347]
NotebookOptionsPosition[     78133,       2319]
NotebookOutlinePosition[     78531,       2335]
CellTagsIndexPosition[     78488,       2332]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[{
 StyleBox["This program evaluates the analytical expressions describing the \
surface uplift above a vertical edge dislocation located a specified distance \
below the surface. Two of these dislocations are then superposed to model \
double ridges on Europa. We assume that Europa\[CloseCurlyQuote]s shell is \
composed of pure water ice. \n\nVariables can be easily changed to consider \
more cases. \n\nThis program is adapted from the solutions presented by \
Savage (1998), please see their work for details.  ", "Subtitle"],
 "\n"
}], "Subtitle",
 CellChangeTimes->{{3.8920429745103283`*^9, 3.892043043456821*^9}, {
  3.892043075719832*^9, 3.892043076237084*^9}, {3.892043136328532*^9, 
  3.892043210672483*^9}, {3.892043244814172*^9, 3.892043261541287*^9}, {
  3.892043292597796*^9, 3.892043294556757*^9}, {3.8920434944948273`*^9, 
  3.892043718068576*^9}, {3.8920437517704372`*^9, 3.892043755133531*^9}, {
  3.892043785922056*^9, 3.892043794144945*^9}, {3.892043869006652*^9, 
  3.892043890210306*^9}, {3.8920465643591948`*^9, 
  3.892046565585256*^9}},ExpressionUUID->"d8ab56f1-2d68-41bb-82cf-\
114c911a01b5"],

Cell[CellGroupData[{

Cell["\<\
First, evaluate the surface displacement caused by an upward facing \
dislocation. 

The ice wedge is:
500 m wide (bv = 5e5 mm)
500 m below surface (c = 500 m = 5e5 mm) 
in a 1 km thick layer of ice (h = 1e6 mm)

Running the cell below will output two plots that show the displacement \
caused by the dislocation in the horizontal and vertical directions, \
respectively. \
\>", "Subsection",
 CellChangeTimes->{
  3.8920433347255907`*^9, {3.892043367709*^9, 3.8920434115553703`*^9}, {
   3.892043948267002*^9, 3.892044009810231*^9}, {3.892044067808873*^9, 
   3.892044154279907*^9}, {3.892044333970284*^9, 3.892044436189612*^9}, {
   3.892046088877108*^9, 
   3.892046102623213*^9}},ExpressionUUID->"06b5dc72-833f-4a81-ab2b-\
0251805ad19e"],

Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{
   RowBox[{
    RowBox[{"bv", "=", "5*^5"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"gam", "=", "0.5"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"kap1", "=", 
     RowBox[{"5", "/", "3"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"kap2", "=", 
     RowBox[{"5", "/", "3"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"c", "=", "5*^5"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"h", "=", "1*^6"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y1", "=", 
     RowBox[{"h", "/", "5"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y2", "=", 
     RowBox[{"2", "*", 
      RowBox[{"h", "/", "5"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y3", "=", 
     RowBox[{"3", "*", 
      RowBox[{"h", "/", "5"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y4", "=", "h"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y5", "=", 
     RowBox[{"3", "*", 
      RowBox[{"h", "/", "2"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y6", "=", 
     RowBox[{"2", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y7", "=", 
     RowBox[{"3", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y8", "=", 
     RowBox[{"5", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y9", "=", 
     RowBox[{"7", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y10", "=", 
     RowBox[{"10", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"sig", "=", "1"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"del", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"gam", "*", "kap1"}], "+", "1"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"kap2", "+", "1"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"eta", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"gam", "+", "kap2"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"kap2", "+", "1"}], ")"}]}]}], ";", 
    RowBox[{"yref", "=", 
     RowBox[{"30", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"k", "=", 
     RowBox[{"c", "/", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"p", "=", 
     RowBox[{"Cosh", "[", 
      RowBox[{"xi", "*", "h"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"q", "=", 
     RowBox[{"Sinh", "[", 
      RowBox[{"xi", "*", "h"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"t", "=", 
     RowBox[{"xi", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"m", "=", 
     RowBox[{"Exp", "[", 
      RowBox[{"2", "*", "k", "*", "t"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"g2", "=", 
     RowBox[{"gam", "*", "g1"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"coef", "=", 
     RowBox[{"2", "*", 
      RowBox[{"(", 
       RowBox[{"del", "+", "eta", "-", "1"}], ")"}], "*", "g1", "*", 
      RowBox[{"bv", "/", 
       RowBox[{"(", 
        RowBox[{"del", "*", "eta", "*", 
         RowBox[{"(", 
          RowBox[{"kap1", "+", "1"}], ")"}]}], ")"}]}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"phi1", "=", 
     RowBox[{"coef", "*", 
      RowBox[{"(", 
       RowBox[{"I", "*", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c", "+", "x"}], ")"}]}], "*", "xi"}], "]"}], 
         "/", "xi"}]}], ")"}], "*", "ca"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"psi1", "=", 
     RowBox[{"coef", "*", 
      RowBox[{"(", 
       RowBox[{"I", "*", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c", "+", "x"}], ")"}]}], "*", "xi"}], "]"}], 
         "/", 
         RowBox[{"xi", "^", "2"}]}]}], ")"}], "*", "cb"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"phi2", "=", 
     RowBox[{"coef", "*", 
      RowBox[{"(", 
       RowBox[{"I", "*", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}]}], "*", "xi"}], "]"}], "/", 
         "xi"}]}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"cc", "*", 
         RowBox[{"Cosh", "[", 
          RowBox[{"xi", "*", "x"}], "]"}]}], "+", "\n", 
        "\[NonBreakingSpace] ", 
        RowBox[{"cd", "*", 
         RowBox[{"Sinh", "[", 
          RowBox[{"xi", "*", "x"}], "]"}]}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"psi2", "=", 
     RowBox[{"coef", "*", 
      RowBox[{"(", 
       RowBox[{"I", "*", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}]}], "*", "xi"}], "]"}], "/", 
         RowBox[{"xi", "^", "2"}]}]}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"ce", "*", 
         RowBox[{"Cosh", "[", 
          RowBox[{"xi", "*", "x"}], "]"}]}], "+", "\n", 
        "\[NonBreakingSpace] ", 
        RowBox[{"cf", "*", 
         RowBox[{"Sinh", "[", 
          RowBox[{"xi", "*", "x"}], "]"}]}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"ux1", "=", 
     RowBox[{
      RowBox[{"Simplify", "[", 
       RowBox[{
        RowBox[{"phi1", "*", 
         RowBox[{"(", 
          RowBox[{"kap1", "+", "1"}], ")"}]}], "-", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"x", "*", "phi1"}], "+", "psi1"}], ",", "x"}], "]"}]}], 
       "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "g1"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uy1", "=", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Simplify", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"x", "*", "phi1"}], "+", "psi1"}], ")"}], "*", "I", "*", 
         "xi"}], "]"}]}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "g1"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"ux2", "=", 
     RowBox[{
      RowBox[{"Simplify", "[", 
       RowBox[{
        RowBox[{"phi2", "*", 
         RowBox[{"(", 
          RowBox[{"kap2", "+", "1"}], ")"}]}], "-", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"x", "*", "phi2"}], "+", "psi2"}], ",", "x"}], "]"}]}], 
       "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "g2"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uy2", "=", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Simplify", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"x", "*", "phi2"}], "+", "psi2"}], ")"}], "*", "I", "*", 
         "xi"}], "]"}]}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "g2"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"equ1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"gam", "*", "kap1", "*", "ca"}], "+", 
       RowBox[{"gam", "*", "cb"}], "-", 
       RowBox[{"kap2", "*", "cc"}], "+", "cf"}], "==", "0"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"equ2", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"gam", "*", "sig", "*", "cb"}], "-", "ce"}], "==", "0"}]}], 
    ";"}], "\n", 
   RowBox[{
    RowBox[{"equ3", "=", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"kap1", "+", "1"}], ")"}], "*", "sig", "*", "ca"}], "+", 
       RowBox[{"2", "*", "sig", "*", "cb"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"kap2", "+", "1"}], ")"}], "*", "cd"}], "-", 
       RowBox[{"2", "*", "ce"}]}], "==", "0"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"equ4", "=", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"kap1", "-", "1"}], ")"}], "*", "ca"}], "+", 
       RowBox[{"2", "*", "cb"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"kap2", "-", "1"}], ")"}], "*", "cc"}], "+", 
       RowBox[{"2", "*", "cf"}]}], "==", "0"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"equ5", "=", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"kap2", "+", "1"}], ")"}], "*", "q"}], "-", 
          RowBox[{"2", "*", "t", "*", "p"}]}], ")"}], "*", "cc"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"kap2", "+", "1"}], ")"}], "*", "p"}], "-", 
          RowBox[{"2", "*", "t", "*", "q"}]}], ")"}], "*", "cd"}], "+", "\n", 
       " ", 
       RowBox[{"2", "*", "p", "*", "ce"}], "-", 
       RowBox[{"2", "*", "q", "*", "cf"}]}], "==", 
      RowBox[{"I", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "*", "del", "*", "eta", "*", "t", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "k"}], ")"}], "*", "m"}], "+", "del", "-", "eta",
          "+", "\n", " ", 
         RowBox[{"2", "*", "del", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "eta"}], ")"}], "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "k"}], ")"}], "*", "t", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", 
            RowBox[{"2", "*", "k", "*", "t"}]}], ")"}]}]}], ")"}]}]}]}], 
    ";"}], "\n", 
   RowBox[{
    RowBox[{"equ6", "=", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"kap2", "-", "1"}], ")"}], "*", "p"}], "-", 
          RowBox[{"2", "*", "t", "*", "q"}]}], ")"}], "*", "cc"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"kap2", "-", "1"}], ")"}], "*", "q"}], "-", 
          RowBox[{"2", "*", "t", "*", "p"}]}], ")"}], "*", "cd"}], "+", "\n", 
       " ", 
       RowBox[{"2", "*", "q", "*", "ce"}], "-", 
       RowBox[{"2", "*", "p", "*", "cf"}]}], "==", 
      RowBox[{
       RowBox[{"-", "I"}], "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "*", "del", "*", "eta", "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1", "-", "k"}], ")"}], "*", "t"}]}], ")"}], "*", "m"}],
          "+", 
         RowBox[{"2", "*", "del", "*", "eta"}], "-", "\n", " ", "del", "-", 
         "eta", "+", 
         RowBox[{"2", "*", "del", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "eta"}], ")"}], "*", 
          RowBox[{"(", 
           RowBox[{"1", "+", "k"}], ")"}], "*", "t"}], "-", 
         RowBox[{"4", "*", "del", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "eta"}], ")"}], "*", "k", "*", 
          RowBox[{"t", "^", "2"}]}]}], ")"}]}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "equ1", ",", "equ2", ",", "equ3", ",", "equ4", ",", "equ5", ",", 
        "equ6"}], "}"}], ",", "\n", " ", 
      RowBox[{"{", 
       RowBox[{"ca", ",", "cb", ",", "cc", ",", "cd", ",", "ce", ",", "cf"}], 
       "}"}]}], "]"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"ux2e", "=", 
     RowBox[{"ux2", " ", "/.", " ", "%"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uy2e", "=", 
     RowBox[{"uy2", " ", "/.", " ", "%%"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"ux2h", "=", 
     RowBox[{"First", "[", 
      RowBox[{"Simplify", "[", 
       RowBox[{"ux2e", " ", "/.", " ", 
        RowBox[{"x", "\[Rule]", 
         RowBox[{"-", "h"}]}]}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uy2h", "=", 
     RowBox[{"First", "[", 
      RowBox[{"Simplify", "[", 
       RowBox[{"uy2e", " ", "/.", " ", 
        RowBox[{"x", "\[Rule]", 
         RowBox[{"-", "h"}]}]}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"daty", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y1", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y1"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y2", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y2"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y3", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y3"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y4", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y4"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y5", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y5"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y6", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y6"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y7", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y7"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y8", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y8"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y9", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y9"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y10", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y10"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"datx", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y1", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y1"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y2", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y2"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y3", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y3"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y4", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y4"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y5", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y5"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y6", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y6"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y7", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y7"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y8", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y8"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y9", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y9"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y10", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y10"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"r1", "=", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"h", "-", "c"}], ")"}], "^", "2"}], "+", 
       RowBox[{"y", "^", "2"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"r2", "=", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "+", 
       RowBox[{"y", "^", "2"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"theta1", "=", 
     RowBox[{"ArcTan", "[", 
      RowBox[{
       RowBox[{"h", "-", "c"}], ",", 
       RowBox[{"-", "y"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"theta2", "=", 
     RowBox[{"ArcTan", "[", 
      RowBox[{
       RowBox[{"h", "+", "c"}], ",", 
       RowBox[{"-", "y"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"a", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", 
        RowBox[{"1", "/", "gam"}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"kap2", "/", "gam"}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"b", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"kap1", "-", 
        RowBox[{"kap2", "/", "gam"}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"kap1", "+", 
        RowBox[{"1", "/", "gam"}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uxyhalf", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"bv", "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", "Pi", 
          RowBox[{"(", 
           RowBox[{"kap2", "+", "1"}], ")"}]}], ")"}]}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "-", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", 
            RowBox[{"Log", "[", "r1", "]"}]}], "-", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "-", "c"}], ")"}], "^", "2"}], "/", 
            RowBox[{"r1", "^", "2"}]}]}], ")"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"b", "+", "a"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "-", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", 
            RowBox[{"Log", "[", "r2", "]"}]}], "-", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "/", 
            RowBox[{"r2", "^", "2"}]}]}], ")"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"b", "-", "a"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "+", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", 
            RowBox[{"Log", "[", "r2", "]"}]}], "-", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "/", 
            RowBox[{"r2", "^", "2"}]}]}], ")"}]}], "+", 
        RowBox[{"2", "*", "c", "*", "a", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "2"}], "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "/", "\n", 
             "\[NonBreakingSpace] ", 
             RowBox[{"r2", "^", "2"}]}]}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"3", "-", "kap2"}], ")"}], "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "/", 
             RowBox[{"r2", "^", "2"}]}]}], "-", 
           RowBox[{"4", "*", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"h", "+", "c"}], ")"}], "^", "3"}], "/", "\n", 
             "\[NonBreakingSpace] ", 
             RowBox[{"r2", "^", "4"}]}]}], "+", 
           RowBox[{"2", "*", "c", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "1"}], "/", 
               RowBox[{"r2", "^", "2"}]}], "+", 
              RowBox[{"2", "*", 
               RowBox[{
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "/", 
                RowBox[{"r2", "^", "4"}]}]}]}], ")"}]}]}], ")"}]}]}], 
       ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uyyhalf", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"bv", "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", "Pi", 
          RowBox[{"(", 
           RowBox[{"kap2", "+", "1"}], ")"}]}], ")"}]}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "+", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", "theta1"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"h", "-", "c"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r1", "^", "2"}]}]}]}], ")"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"b", "+", "a"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "+", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", "theta2"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "2"}]}]}]}], ")"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"b", "-", "a"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "-", "\n", "\[NonBreakingSpace] ", "1"}], ")"}], 
            "*", "theta2"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "2"}]}]}]}], ")"}]}], "+", 
        RowBox[{"2", "*", "a", "*", "c", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "*", 
            RowBox[{"y", "/", "\n", "\[NonBreakingSpace] ", 
             RowBox[{"r2", "^", "2"}]}]}], "-", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"1", "+", "kap2"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "2"}]}]}], "+", 
           RowBox[{"4", "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "4"}]}]}], "\n", "\[NonBreakingSpace] ", "-", 
           
           RowBox[{"4", "*", "c", "*", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "4"}]}]}]}], ")"}]}]}], ")"}]}]}], ";"}], 
   "\n", 
   RowBox[{
    RowBox[{"fidxy", "=", 
     RowBox[{"N", "[", 
      RowBox[{"uxyhalf", " ", "/.", " ", 
       RowBox[{"y", "\[Rule]", "0"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"fidyy", "=", 
     RowBox[{"N", "[", 
      RowBox[{"uyyhalf", " ", "/.", " ", 
       RowBox[{"y", "\[Rule]", "0"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"datxyhalf", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y1"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y2"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y3"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y4"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y5"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y6"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y7"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y8"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y9"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y10"}]}], "]"}], "-", "fidxy"}]}], 
        "}"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"datyyhalf", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y1"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y2"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y3"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y4"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y5"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y6"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y7"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y8"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y9"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y10"}]}], "]"}], "-", "fidyy"}]}], 
        "}"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{"datxldh", "=", 
    RowBox[{"datx", "+", "datxyhalf"}]}], "\n", 
   RowBox[{"datyldh", "=", 
    RowBox[{"daty", "+", "datyyhalf"}]}], "\n", 
   RowBox[{"uxfin", "=", 
    RowBox[{"ListPlot", "[", 
     RowBox[{
      RowBox[{"datx", "+", "datxyhalf"}], ",", 
      RowBox[{"PlotStyle", "\[Rule]", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"PointSize", "[", "0.015", "]"}], "}"}]}], ",", 
      RowBox[{"PlotRange", "\[Rule]", "All"}], ",", " ", 
      RowBox[{
      "PlotLabel", "\[Rule]", "\"\<Horizontal displacement (Ux)\>\""}], ",", 
      " ", 
      RowBox[{"AxesLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{"mm", ",", " ", "mm"}], "}"}]}]}], "]"}]}], "\n", 
   RowBox[{"uyfin", "=", 
    RowBox[{"ListPlot", "[", 
     RowBox[{
      RowBox[{"daty", "+", "datyyhalf"}], ",", 
      RowBox[{"PlotStyle", "\[Rule]", "\n", 
       RowBox[{"{", 
        RowBox[{"PointSize", "[", "0.02", "]"}], "}"}]}], ",", 
      RowBox[{"PlotRange", "\[Rule]", "All"}], ",", " ", 
      RowBox[{"PlotLabel", "\[Rule]", "\"\<Vertical displacement (Uy)\>\""}], 
      ",", " ", 
      RowBox[{"AxesLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{"mm", ",", " ", "mm"}], "}"}]}]}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{
  3.8920429745103283`*^9, {3.8920461073753567`*^9, 3.892046110196101*^9}, {
   3.892046459071898*^9, 
   3.892046473087296*^9}},ExpressionUUID->"a66c8ad1-09b5-47d5-a617-\
477085516254"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Now evaluate the displacement caused by a downward facing dislocation. 

The ice wedge is:
500 m wide but facing the opposite direction as before (bv = -5e5 mm)
100 m below surface (c = 900 m = 9e5 mm) 
in a 1 km thick layer of ice (h = 1e6 mm)

Running the cell below will output two plots that show the displacement \
caused by the dislocation in the horizontal and vertical directions, \
respectively. \
\>", "Subsection",
 CellChangeTimes->{{3.892044584362048*^9, 3.8920447294509993`*^9}, {
  3.892045884563035*^9, 
  3.892045937532889*^9}},ExpressionUUID->"98a3ed96-b565-48ac-a0c6-\
75f3c6d82de6"],

Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{
   RowBox[{
    RowBox[{"bv", "=", 
     RowBox[{"-", "5*^5"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"gam", "=", "0.5"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"kap1", "=", 
     RowBox[{"5", "/", "3"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"kap2", "=", 
     RowBox[{"5", "/", "3"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"c", "=", "9.0*^5"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"h", "=", "1*^6"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y1", "=", 
     RowBox[{"h", "/", "5"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y2", "=", 
     RowBox[{"2", "*", 
      RowBox[{"h", "/", "5"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y3", "=", 
     RowBox[{"3", "*", 
      RowBox[{"h", "/", "5"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y4", "=", "h"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y5", "=", 
     RowBox[{"3", "*", 
      RowBox[{"h", "/", "2"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y6", "=", 
     RowBox[{"2", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y7", "=", 
     RowBox[{"3", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y8", "=", 
     RowBox[{"5", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y9", "=", 
     RowBox[{"7", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"y10", "=", 
     RowBox[{"10", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"sig", "=", "1"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"del", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"gam", "*", "kap1"}], "+", "1"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"kap2", "+", "1"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"eta", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"gam", "+", "kap2"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"kap2", "+", "1"}], ")"}]}]}], ";", 
    RowBox[{"yref", "=", 
     RowBox[{"30", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"k", "=", 
     RowBox[{"c", "/", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"p", "=", 
     RowBox[{"Cosh", "[", 
      RowBox[{"xi", "*", "h"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"q", "=", 
     RowBox[{"Sinh", "[", 
      RowBox[{"xi", "*", "h"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"t", "=", 
     RowBox[{"xi", "*", "h"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"m", "=", 
     RowBox[{"Exp", "[", 
      RowBox[{"2", "*", "k", "*", "t"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"g2", "=", 
     RowBox[{"gam", "*", "g1"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"coef", "=", 
     RowBox[{"2", "*", 
      RowBox[{"(", 
       RowBox[{"del", "+", "eta", "-", "1"}], ")"}], "*", "g1", "*", 
      RowBox[{"bv", "/", 
       RowBox[{"(", 
        RowBox[{"del", "*", "eta", "*", 
         RowBox[{"(", 
          RowBox[{"kap1", "+", "1"}], ")"}]}], ")"}]}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"phi1", "=", 
     RowBox[{"coef", "*", 
      RowBox[{"(", 
       RowBox[{"I", "*", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c", "+", "x"}], ")"}]}], "*", "xi"}], "]"}], 
         "/", "xi"}]}], ")"}], "*", "ca"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"psi1", "=", 
     RowBox[{"coef", "*", 
      RowBox[{"(", 
       RowBox[{"I", "*", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c", "+", "x"}], ")"}]}], "*", "xi"}], "]"}], 
         "/", 
         RowBox[{"xi", "^", "2"}]}]}], ")"}], "*", "cb"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"phi2", "=", 
     RowBox[{"coef", "*", 
      RowBox[{"(", 
       RowBox[{"I", "*", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}]}], "*", "xi"}], "]"}], "/", 
         "xi"}]}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"cc", "*", 
         RowBox[{"Cosh", "[", 
          RowBox[{"xi", "*", "x"}], "]"}]}], "+", "\n", 
        "\[NonBreakingSpace] ", 
        RowBox[{"cd", "*", 
         RowBox[{"Sinh", "[", 
          RowBox[{"xi", "*", "x"}], "]"}]}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"psi2", "=", 
     RowBox[{"coef", "*", 
      RowBox[{"(", 
       RowBox[{"I", "*", 
        RowBox[{
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}]}], "*", "xi"}], "]"}], "/", 
         RowBox[{"xi", "^", "2"}]}]}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"ce", "*", 
         RowBox[{"Cosh", "[", 
          RowBox[{"xi", "*", "x"}], "]"}]}], "+", "\n", 
        "\[NonBreakingSpace] ", 
        RowBox[{"cf", "*", 
         RowBox[{"Sinh", "[", 
          RowBox[{"xi", "*", "x"}], "]"}]}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"ux1", "=", 
     RowBox[{
      RowBox[{"Simplify", "[", 
       RowBox[{
        RowBox[{"phi1", "*", 
         RowBox[{"(", 
          RowBox[{"kap1", "+", "1"}], ")"}]}], "-", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"x", "*", "phi1"}], "+", "psi1"}], ",", "x"}], "]"}]}], 
       "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "g1"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uy1", "=", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Simplify", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"x", "*", "phi1"}], "+", "psi1"}], ")"}], "*", "I", "*", 
         "xi"}], "]"}]}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "g1"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"ux2", "=", 
     RowBox[{
      RowBox[{"Simplify", "[", 
       RowBox[{
        RowBox[{"phi2", "*", 
         RowBox[{"(", 
          RowBox[{"kap2", "+", "1"}], ")"}]}], "-", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"x", "*", "phi2"}], "+", "psi2"}], ",", "x"}], "]"}]}], 
       "]"}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "g2"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uy2", "=", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Simplify", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"x", "*", "phi2"}], "+", "psi2"}], ")"}], "*", "I", "*", 
         "xi"}], "]"}]}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "g2"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"equ1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"gam", "*", "kap1", "*", "ca"}], "+", 
       RowBox[{"gam", "*", "cb"}], "-", 
       RowBox[{"kap2", "*", "cc"}], "+", "cf"}], "==", "0"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"equ2", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"gam", "*", "sig", "*", "cb"}], "-", "ce"}], "==", "0"}]}], 
    ";"}], "\n", 
   RowBox[{
    RowBox[{"equ3", "=", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"kap1", "+", "1"}], ")"}], "*", "sig", "*", "ca"}], "+", 
       RowBox[{"2", "*", "sig", "*", "cb"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"kap2", "+", "1"}], ")"}], "*", "cd"}], "-", 
       RowBox[{"2", "*", "ce"}]}], "==", "0"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"equ4", "=", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"kap1", "-", "1"}], ")"}], "*", "ca"}], "+", 
       RowBox[{"2", "*", "cb"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"kap2", "-", "1"}], ")"}], "*", "cc"}], "+", 
       RowBox[{"2", "*", "cf"}]}], "==", "0"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"equ5", "=", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"kap2", "+", "1"}], ")"}], "*", "q"}], "-", 
          RowBox[{"2", "*", "t", "*", "p"}]}], ")"}], "*", "cc"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"kap2", "+", "1"}], ")"}], "*", "p"}], "-", 
          RowBox[{"2", "*", "t", "*", "q"}]}], ")"}], "*", "cd"}], "+", "\n", 
       " ", 
       RowBox[{"2", "*", "p", "*", "ce"}], "-", 
       RowBox[{"2", "*", "q", "*", "cf"}]}], "==", 
      RowBox[{"I", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "*", "del", "*", "eta", "*", "t", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "k"}], ")"}], "*", "m"}], "+", "del", "-", "eta",
          "+", "\n", " ", 
         RowBox[{"2", "*", "del", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "eta"}], ")"}], "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "k"}], ")"}], "*", "t", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", 
            RowBox[{"2", "*", "k", "*", "t"}]}], ")"}]}]}], ")"}]}]}]}], 
    ";"}], "\n", 
   RowBox[{
    RowBox[{"equ6", "=", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"kap2", "-", "1"}], ")"}], "*", "p"}], "-", 
          RowBox[{"2", "*", "t", "*", "q"}]}], ")"}], "*", "cc"}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"kap2", "-", "1"}], ")"}], "*", "q"}], "-", 
          RowBox[{"2", "*", "t", "*", "p"}]}], ")"}], "*", "cd"}], "+", "\n", 
       " ", 
       RowBox[{"2", "*", "q", "*", "ce"}], "-", 
       RowBox[{"2", "*", "p", "*", "cf"}]}], "==", 
      RowBox[{
       RowBox[{"-", "I"}], "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "*", "del", "*", "eta", "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1", "-", "k"}], ")"}], "*", "t"}]}], ")"}], "*", "m"}],
          "+", 
         RowBox[{"2", "*", "del", "*", "eta"}], "-", "\n", " ", "del", "-", 
         "eta", "+", 
         RowBox[{"2", "*", "del", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "eta"}], ")"}], "*", 
          RowBox[{"(", 
           RowBox[{"1", "+", "k"}], ")"}], "*", "t"}], "-", 
         RowBox[{"4", "*", "del", "*", 
          RowBox[{"(", 
           RowBox[{"1", "-", "eta"}], ")"}], "*", "k", "*", 
          RowBox[{"t", "^", "2"}]}]}], ")"}]}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "equ1", ",", "equ2", ",", "equ3", ",", "equ4", ",", "equ5", ",", 
        "equ6"}], "}"}], ",", "\n", " ", 
      RowBox[{"{", 
       RowBox[{"ca", ",", "cb", ",", "cc", ",", "cd", ",", "ce", ",", "cf"}], 
       "}"}]}], "]"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"ux2e", "=", 
     RowBox[{"ux2", " ", "/.", " ", "%"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uy2e", "=", 
     RowBox[{"uy2", " ", "/.", " ", "%%"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"ux2h", "=", 
     RowBox[{"First", "[", 
      RowBox[{"Simplify", "[", 
       RowBox[{"ux2e", " ", "/.", " ", 
        RowBox[{"x", "\[Rule]", 
         RowBox[{"-", "h"}]}]}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uy2h", "=", 
     RowBox[{"First", "[", 
      RowBox[{"Simplify", "[", 
       RowBox[{"uy2e", " ", "/.", " ", 
        RowBox[{"x", "\[Rule]", 
         RowBox[{"-", "h"}]}]}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"daty", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y1", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y1"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y2", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y2"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y3", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y3"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y4", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y4"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y5", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y5"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y6", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y6"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y7", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y7"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y8", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y8"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y9", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y9"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y10", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"I", "/", "Pi"}], "*", "uy2h", "*", 
            RowBox[{"Sin", "[", 
             RowBox[{"xi", "*", "y10"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"datx", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y1", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y1"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y2", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y2"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y3", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y3"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y4", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y4"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y5", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y5"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y6", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y6"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y7", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y7"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y8", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y8"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y9", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y9"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}], ",", "\n", 
       RowBox[{"{", 
        RowBox[{"y10", ",", 
         RowBox[{"NIntegrate", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"ux2h", "/", "Pi"}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Cos", "[", 
               RowBox[{"xi", "*", "y10"}], "]"}], "-", "1"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"xi", ",", "0", ",", "1", ",", "Infinity"}], "}"}]}], 
          "]"}]}], "}"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"r1", "=", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"h", "-", "c"}], ")"}], "^", "2"}], "+", 
       RowBox[{"y", "^", "2"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"r2", "=", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "+", 
       RowBox[{"y", "^", "2"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"theta1", "=", 
     RowBox[{"ArcTan", "[", 
      RowBox[{
       RowBox[{"h", "-", "c"}], ",", 
       RowBox[{"-", "y"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"theta2", "=", 
     RowBox[{"ArcTan", "[", 
      RowBox[{
       RowBox[{"h", "+", "c"}], ",", 
       RowBox[{"-", "y"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"a", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", 
        RowBox[{"1", "/", "gam"}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"kap2", "/", "gam"}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"b", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"kap1", "-", 
        RowBox[{"kap2", "/", "gam"}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"kap1", "+", 
        RowBox[{"1", "/", "gam"}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uxyhalf", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"bv", "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", "Pi", 
          RowBox[{"(", 
           RowBox[{"kap2", "+", "1"}], ")"}]}], ")"}]}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "-", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", 
            RowBox[{"Log", "[", "r1", "]"}]}], "-", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "-", "c"}], ")"}], "^", "2"}], "/", 
            RowBox[{"r1", "^", "2"}]}]}], ")"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"b", "+", "a"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "-", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", 
            RowBox[{"Log", "[", "r2", "]"}]}], "-", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "/", 
            RowBox[{"r2", "^", "2"}]}]}], ")"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"b", "-", "a"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "+", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", 
            RowBox[{"Log", "[", "r2", "]"}]}], "-", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "/", 
            RowBox[{"r2", "^", "2"}]}]}], ")"}]}], "+", 
        RowBox[{"2", "*", "c", "*", "a", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "2"}], "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "/", "\n", 
             "\[NonBreakingSpace] ", 
             RowBox[{"r2", "^", "2"}]}]}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"3", "-", "kap2"}], ")"}], "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "/", 
             RowBox[{"r2", "^", "2"}]}]}], "-", 
           RowBox[{"4", "*", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"h", "+", "c"}], ")"}], "^", "3"}], "/", "\n", 
             "\[NonBreakingSpace] ", 
             RowBox[{"r2", "^", "4"}]}]}], "+", 
           RowBox[{"2", "*", "c", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "1"}], "/", 
               RowBox[{"r2", "^", "2"}]}], "+", 
              RowBox[{"2", "*", 
               RowBox[{
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "/", 
                RowBox[{"r2", "^", "4"}]}]}]}], ")"}]}]}], ")"}]}]}], 
       ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"uyyhalf", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"bv", "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", "Pi", 
          RowBox[{"(", 
           RowBox[{"kap2", "+", "1"}], ")"}]}], ")"}]}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "+", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", "theta1"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"h", "-", "c"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r1", "^", "2"}]}]}]}], ")"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"b", "+", "a"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "+", "1"}], ")"}], "*", "\n", 
            "\[NonBreakingSpace] ", "theta2"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "2"}]}]}]}], ")"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"b", "-", "a"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"0.5", "*", 
            RowBox[{"(", 
             RowBox[{"kap2", "-", "\n", "\[NonBreakingSpace] ", "1"}], ")"}], 
            "*", "theta2"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "2"}]}]}]}], ")"}]}], "+", 
        RowBox[{"2", "*", "a", "*", "c", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "*", 
            RowBox[{"y", "/", "\n", "\[NonBreakingSpace] ", 
             RowBox[{"r2", "^", "2"}]}]}], "-", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"1", "+", "kap2"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "2"}]}]}], "+", 
           RowBox[{"4", "*", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"h", "+", "c"}], ")"}], "^", "2"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "4"}]}]}], "\n", "\[NonBreakingSpace] ", "-", 
           
           RowBox[{"4", "*", "c", "*", 
            RowBox[{"(", 
             RowBox[{"h", "+", "c"}], ")"}], "*", 
            RowBox[{"y", "/", 
             RowBox[{"r2", "^", "4"}]}]}]}], ")"}]}]}], ")"}]}]}], ";"}], 
   "\n", 
   RowBox[{
    RowBox[{"fidxy", "=", 
     RowBox[{"N", "[", 
      RowBox[{"uxyhalf", " ", "/.", " ", 
       RowBox[{"y", "\[Rule]", "0"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"fidyy", "=", 
     RowBox[{"N", "[", 
      RowBox[{"uyyhalf", " ", "/.", " ", 
       RowBox[{"y", "\[Rule]", "0"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"datxyhalf", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y1"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y2"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y3"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y4"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y5"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y6"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y7"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y8"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y9"}]}], "]"}], "-", "fidxy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uxyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y10"}]}], "]"}], "-", "fidxy"}]}], 
        "}"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"datyyhalf", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y1"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y2"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y3"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y4"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y5"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y6"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y7"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y8"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y9"}]}], "]"}], "-", "fidyy"}]}], "}"}], 
       ",", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"uyyhalf", " ", "/.", " ", 
            RowBox[{"y", "\[Rule]", "y10"}]}], "]"}], "-", "fidyy"}]}], 
        "}"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{"datxldh", "=", 
    RowBox[{"datx", "+", "datxyhalf"}]}], "\n", 
   RowBox[{"datyldh", "=", 
    RowBox[{"daty", "+", "datyyhalf"}]}], "\n", 
   RowBox[{"uxfin", "=", 
    RowBox[{"ListPlot", "[", 
     RowBox[{
      RowBox[{"datx", "+", "datxyhalf"}], ",", 
      RowBox[{"PlotStyle", "\[Rule]", "\n", " ", 
       RowBox[{"{", 
        RowBox[{"PointSize", "[", "0.015", "]"}], "}"}]}], ",", 
      RowBox[{"PlotRange", "\[Rule]", "All"}], ",", " ", 
      RowBox[{
      "PlotLabel", "\[Rule]", "\"\<Horizontal displacement (Ux)\>\""}], ",", 
      " ", 
      RowBox[{"AxesLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{"mm", ",", " ", "mm"}], "}"}]}]}], "]"}]}], "\n", 
   RowBox[{"uyfin", "=", 
    RowBox[{"ListPlot", "[", 
     RowBox[{
      RowBox[{"daty", "+", "datyyhalf"}], ",", 
      RowBox[{"PlotStyle", "\[Rule]", "\n", 
       RowBox[{"{", 
        RowBox[{"PointSize", "[", "0.02", "]"}], "}"}]}], ",", 
      RowBox[{"PlotRange", "\[Rule]", "All"}], ",", " ", 
      RowBox[{"PlotLabel", "\[Rule]", "\"\<Vertical displacement (Uy)\>\""}], 
      ",", " ", 
      RowBox[{"AxesLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{"mm", ",", " ", "mm"}], "}"}]}]}], "]"}]}], "\n"}]}]], "Input",\

 CellChangeTimes->{{3.892044592103404*^9, 3.892044595904217*^9}, {
  3.892045889030342*^9, 3.892045889197857*^9}, {3.892045955861003*^9, 
  3.89204596093294*^9}, {3.8920463083934813`*^9, 3.892046396258423*^9}},
 CellLabel->"In[208]:=",
 CloudRendering -> 
 "InteractiveBoxes",ExpressionUUID->"0e23c80f-f51c-4391-b1e9-616f548b6d8d"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Ridges are modeled as the superposition of the vertical displacements from \
both dislocations. 

These opposing dislocations are at depths of 100 m and 500 m and have bv= \
|500 m|\
\>", "Subsection",
 CellChangeTimes->{{3.892044828690771*^9, 3.8920448707598667`*^9}, {
  3.892044904967729*^9, 3.892044933146401*^9}, {3.892044976894628*^9, 
  3.892044997793515*^9}, {3.8920450287296057`*^9, 3.89204503047285*^9}, {
  3.892045333215251*^9, 
  3.8920454189898663`*^9}},ExpressionUUID->"df9df018-4221-46c0-8617-\
a86ffd79a13f"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"top", " ", "=", " ", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"200000", ",", 
        RowBox[{"66384.7", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"400000", ",", 
        RowBox[{"95656.2", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"600000", ",", 
        RowBox[{"108222.", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1000000", ",", 
        RowBox[{"120255.", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1500000", ",", 
        RowBox[{"127778.", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"2000000", ",", 
        RowBox[{"132562.", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"3000000", ",", 
        RowBox[{"139145.", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"5000000", ",", 
        RowBox[{"147015.", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"7000000", ",", 
        RowBox[{"151366.", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"10000000", ",", 
        RowBox[{"155006.", "\[VeryThinSpace]", "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}]}], "}"}]}], ";"}], 
  "\n"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"bottom", " ", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"200000", ",", 
        RowBox[{
         RowBox[{"-", "12389.2"}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"400000", ",", 
        RowBox[{
         RowBox[{"-", "29761.1"}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"600000", ",", 
        RowBox[{
         RowBox[{"-", "47725.3"}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1000000", ",", 
        RowBox[{
         RowBox[{"-", "75098."}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1500000", ",", 
        RowBox[{
         RowBox[{"-", "95802.6"}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"2000000", ",", 
        RowBox[{
         RowBox[{"-", "108719."}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"3000000", ",", 
        RowBox[{
         RowBox[{"-", "124151."}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"5000000", ",", 
        RowBox[{
         RowBox[{"-", "138842."}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"7000000", ",", 
        RowBox[{
         RowBox[{"-", "145795."}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"10000000", ",", 
        RowBox[{
         RowBox[{"-", "151225."}], "+", 
         RowBox[{"0.", " ", "\[ImaginaryI]"}]}]}], "}"}]}], "}"}]}], ";"}], 
  "\n"}], "\n", 
 RowBox[{
  RowBox[{"bothx", " ", "=", 
   RowBox[{"1*^-6", "*", 
    RowBox[{"top", "[", 
     RowBox[{"[", 
      RowBox[{"All", ",", "1"}], "]"}], "]"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"bothy", " ", "=", 
   RowBox[{"1*^-3", "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"top", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "2"}], "]"}], "]"}], "+", 
      RowBox[{"bottom", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ")"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"both", " ", "=", 
    RowBox[{"Transpose", "[", 
     RowBox[{"{", 
      RowBox[{"bothx", ",", " ", "bothy"}], "}"}], "]"}]}], ";"}], 
  "\n"}], "\n", 
 RowBox[{"fin", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"both", ",", 
    RowBox[{"PlotStyle", "\[Rule]", "\n", " ", 
     RowBox[{"{", 
      RowBox[{"PointSize", "[", "0.015", "]"}], "}"}]}], ",", 
    RowBox[{"PlotRange", "\[Rule]", "All"}], ",", " ", 
    RowBox[{
    "PlotLabel", "\[Rule]", 
     "\"\<Vertical displacement (Uy) from both dislocations\>\""}], ",", " ", 
    
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"km", ",", " ", "meters"}], "}"}]}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.8920453097401667`*^9, 3.892045310799431*^9}, {
  3.892045466912126*^9, 3.892045469622547*^9}, {3.892045523934177*^9, 
  3.892045538871311*^9}},
 CellLabel->"In[76]:=",ExpressionUUID->"316ad221-3500-4c12-8332-46270007d400"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1268, 881},
WindowMargins->{{180, Automatic}, {Automatic, 22}},
FrontEndVersion->"13.2 for Mac OS X x86 (64-bit) (January 30, 2023)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"d1f77760-c022-459a-a59a-080b74204cad"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 1133, 18, 277, "Subtitle",ExpressionUUID->"d8ab56f1-2d68-41bb-82cf-114c911a01b5"],
Cell[CellGroupData[{
Cell[1738, 44, 751, 19, 270, "Subsection",ExpressionUUID->"06b5dc72-833f-4a81-ab2b-0251805ad19e"],
Cell[2492, 65, 34489, 1029, 2383, "Input",ExpressionUUID->"a66c8ad1-09b5-47d5-a617-477085516254"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37018, 1099, 612, 15, 270, "Subsection",ExpressionUUID->"98a3ed96-b565-48ac-a0c6-75f3c6d82de6"],
Cell[37633, 1116, 34646, 1032, 2404, "Input",ExpressionUUID->"0e23c80f-f51c-4391-b1e9-616f548b6d8d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[72316, 2153, 535, 12, 108, "Subsection",ExpressionUUID->"df9df018-4221-46c0-8617-a86ffd79a13f"],
Cell[72854, 2167, 5251, 148, 304, "Input",ExpressionUUID->"316ad221-3500-4c12-8332-46270007d400"]
}, Open  ]]
}, Open  ]]
}
]
*)

